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ABSTRACT 


'The  present  paper  explores  a  method  to  locate  fracture  dam¬ 
age  in  a  structure,  and  to  estimate  the  intensity  of  the  damage. 
It  is  noted  that  when  fracture  damage  occurs  in  a  structure,  the 
modal  frquencies  of  the  damaged  structure  will  shift.  The  per¬ 
turbations  in  individual  modal  frequencies  are  governed  uniquely 
by  the  intensity  of  the  damage  and  its  location  in  the  struc¬ 
ture.  Consequently,  when  the  frequency  spectrum  of  an  unflawed 
structure  is  available,  field  measurement  of  the  structural 
response  after  a  major  event  could  lead  to  an  estimate  of  inten¬ 
sity  and  location  of  its  damage.  - 

The  existence  of  a  crack  renders  an  overall  softening  effect 
in  the  compliance  of  the  structure.  The  reduction  in  stiffness 
can  be  quantified  by  the  location  and  the  intensity  of  the 
crack.  With  such  a  model,  the  diagnosis  of  the  fracture  damage 
is  fundamentally  a  structural  identification*  The  effect  of 
fracture  at  a  section  may  be  represented  by  a  Spring-loaded  hinge 
of  a  certain  torsional  spring  coefficient  k.  Such  representation 
is  based  on  the  criteria  of  displacement  continuity,  slope  dis¬ 
continuity  at  the  cracked  section  and  increased  compliance.  It 
is  demonstrated  in  the  report  that  the  spring  coefficient  is  con¬ 
stant  for  cross  cracks.  For  one-sided  or  nonsymmetrical  two- 
sided  cracks,  the  "bilinear  spring"  may  be  approximated  by  a 
linearized  "fracture  spring  hinge".  It  is  also  demonstrated  that 
the  spring  "constant"  is  indeed  invariant  to  the  location  of  the 
fracture.  The  spring  constant  is  computed  by  the  strain  energy 
release  on  the  crack  surface.  The  shifted  modal  frequencies  are 
readily  the  modal  frequencies  of  the  same  structure  with  the 
presence  of  such  a  "fracture  hinge". 

Simple  beams,  l-frames,  and  U-frames  have  been  chosen  as  the 
basic  simple  structures.  For  illustrative  purposes,  graphs  are 
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constructed  in  damage-versus-location  space,  for  a  specific  fre¬ 
quency  variation  which  results  from  a  crack  at  various  locations 
on  the  beam  with  corresponding  intensity  of  damage.  A  composite 
of  such  graphs  for  several  modal  frequency  variations  thus  yields 
the  point  of  specific  location  and  the  intensity  of  the  fracture 
damage.  Analytically,  the  damage  function  has  been  programmed  to 
solve  for  the  fracture  damage  characteristics  (y,  e) ,  signifying 
the  intensity  of  the  damage  $nd  its  location.  When  there  are 
inadequate  measurements  of  modal  frequency  variations,  the  frac¬ 
ture  damage  characteristics  (y,  e)  will  either  be  multiple-valued 
or  be  some  continuous  function  defined  over  the  structure.  Some 
probability  distribution  function  must  be  determined  or  assumed, 
either  from  statistics  or  load  distribution  or  a  combination 
thereof,  to  give  a  probability  measure  to  the  fracture  damage 
occurrence.  A  multiple-valued  case  has  been  illustrated  in  the 
report. 
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NOMENCLATURE 


a  crack  length 

A  one  side  of  crack  surface  area 

constant  coefficients;  specific  locations  on  the 
structure 

b  half  depth  of  a  flexural  member 

e  crack  location  (=  L*/L) 

E  modulus  of  elasticity 

f( y)  dimensionless  fracture  intensity  factor 

g(8)  frequency  equation 

fracture  toughness 
Gk(R-)  damage  function 

I  area  moment  of  inertia 

Kj  stress  intensity  factor  for  Mode  I  (open  mode)  crack 

L  characteristic  length  of  a  beam 

Ui  length  of  ith  members  in  a  structure 

li  length  ratio  =  L/L^ 

M(x)  resisting  bending  moment  (at  section  x)  *  Ely" (x) 

n^  stiffness  ratio  *  (EI)i/(EI) 

Ri  relative  frequency  variations  of  ith  mode  =  Au^/u^ 

U(P)  submatrix  defined  at  the  location  P 

AU  strain  energy  increment 

V(x)  resisting  cross-shear  (at  section  x)  *  -EIy'"(x) 

x  axial  position 


3 


f 


modal  shape  function  of  ith  beam 
beam  width 

lateral  deflection  of  ith  beam 
characteristic  ratio  =  ^/e 
characteristic  value  s  =  ui  pL  /(El) 
fracture  damage  =  a/b 

torsional  spring  constant  of  the  fracture  hinge 

Poisson's  ratio 

lineal  mass  density 

sensitivity  nisnber 

ith  modal  frequency 

ith  modal  frequency  shift 

dimensionless  axial  position  *  x/L 


1.0  INTRODUCTION 

The  problem  of  fracture  diagnosis  arises  from  the  need  to 
assure  integrity  of  structures  within  the  means  of  nondestructive 
testing.  The  present  report  summarizes  the  modal  method  to 
detect  and  to  identify  a  fracture  damage  in  a  simple  structure. 
The  simple  structures  are  chosen  in  the  initial  effort  essen¬ 
tially  for  their  analytical  feasibility  both  in  computations  and 
verifications.  The  very  concept  shall  later  be  extended  to  more 
complex  structures  and  possibly  multiple  fractures. 

It  is  well  known  and  was  studied  experimentally  and  analyti¬ 
cally  that  a  structural  member  loses  its  load-carrying  capacity 
when  a  crack  is  developed  in  the  member  [1-3]. 1  Structures, 
which  are  designed  to  withstand  unusual  severe  loadings  such  as 
earthquakes  to  buildings,  sudden  pressurization  in  reactor  ves¬ 
sels,  or  blast  loads  to  aerospace  structures,  may  develop  cracks 
either  through  aging  or  as  a  result  of  an  event  of  excitation. 

The  damaged  structure,  if  not  detected  and  corrected,  may  not  be 
able  to  withstand  the  severe  loadings  according  to  design.  ASME 
Code  [4]  has  identified  crack  inspection  methods  by  visual,  sur¬ 
face,  and  volumetric  methods.  The  cracks  may  escape  the  visual 
and  surface  examinations  because  of  sheltering  by  coating  or 
other  structural  components,  or  simply  due  to  oversight.  Radio- 
graphic  volumetric  method  has  been  successfully  employed  for  some 
structures.  It  is  however  limited  by  its  expense,  by  the  range 
limitation,  and  by  the  complexity  of  the  structure.  In  searching 
for  other  means,  it  is  proposed  in  the  present  report  that  frac¬ 
ture,  in  general,  changes  the  dynamic  characteristics  of  the 
structure.  Furthermore,  for  small  oscillations,  the  dynamic 
parameters  are  assumed  to  be  delineative  in  the  damping  and  the 
stiffness  properties.  This  report  shall  study  only  the  change  in 
the  stiffness  parameter,  and  its  relationship  to  fracture 
characteristics  of  the  structure.  The  change  is  in  turn  to  be 
measured  by  the  change  in  modal  frequencies. 

Numbers  in  brackets  denote  the  references. 
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Recently,  Chondros  and  Dimarogonas  [5]  experimentally  deter¬ 
mined  the  change  of  modal  frequencies  of  a  cantilever  beam,  which 
was  caused  by  crack  at  the  welded  built-in  end.  The  changes  in 
frequencies  are  related  analytically  then  to  those  of  a  canti¬ 
lever  beam  which  is  supported  by  a  spring-loaded  hinge  at  the 
built-in  end.  In  the  present  report,  the  basic  concept  will  be 
expanded  further  to  establish  a  theoretical  base  of  fracture 
diagnosis  by  measurement  of  frequency  variations.  The  method 
will  then  be  demonstrated  applicable  to  other  configurations  of 
simple  structures  with  the  crack  at  any  possible  location. 
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2.0  THEORY 

2.1  Damage  Function 

The  fractured  structure  is  to  be  characterized  by  the  frac¬ 
ture  intensity  (y)  and  the  location  (e).  The  fracture  intensity 
is  essentially  defined  through  the  geometry  of  the  crack.  If 
there  are  N  cracks,  the  fracture  damage  of  the  structure  is  hence 
defined  by  a  set  of  k  =  1,2,...,N.  Because  of  the 

presence  of  the  N  cracks,  the.  modal  frequencies  {w-j}  will  be 
changed,  usually  reduced,  by  the  amount  { } .  The  set  of 
relative  frequency  variations  {R-j}  is  defined  as 

Ri  =  Aa>.j  / .  (1) 

The  relative  frequency  variation  can  be  functionally  expressed  as 
Ri  =  *VYk’  ek  ’  structura^  parameters),  (2) 

when  it  is  assumed  that  the  frequency  variations  can  be  measured 
immediately  before  and  after  the  damage-causing  event.  The 
structural  parameters  may  be  considered  as  constants.  A  simpli¬ 
fied  form  of  Equation  (2)  is  hence 

Ri  *  MYk’  ek^ '  ^ 

If  we  exclude  the  case  of  repeated  modal  frequencies,  the  inverse 
of  Equation  (3) 2  exists  as 

(y.  e)k  -  Gk(R.),  (4) 

where  every  G|<  represents  a  pair  of  functions  for  the  pair  of 
fracture  characteristics  (y,  e)^. 


Numbers  in  parentheses  refer  to  equations. 
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With  known  and  sufficient  measurements  of  the  frequency 
variations,  R-j ,  the  damage  characteristics  can  be  determined 
with  given  functional  expressions  {G^fR-j)},  which  are  there¬ 
fore  referred  to  as  the  damage  functions.  The  damage  functions 
can  be  established  experimentally,  as  done  by  Chondros  and 
Dimarogonas  [5]  for  a  fixed  location  e  =  0.  They  may  also  be 
determined  analytically;  the  present  report  addresses  this  prob¬ 
lem  in  the  sequel.  Theoretically,  given  adequate  measurements  of 
the  frequency  variations  the  fracture  damage  (y.  e)k  can  be 
determined  uniquely.  But  uncertainty  rises  through  many  factors, 
such  as  the  number  of  measurable  modal  frequencies,  the  possible 
cross-over  of  frequency-shifts.  Even  the  damage  functions  them¬ 
selves  may  be  probabilistic.  As  a  general  illustration,  let 
there  be  inadequate  M  measurements  of  frequency  variations  for  a 
known  number,  N,  of  cracks,  such  that  M  <  2N  +  1.  The  damage 
function  (4)  will  not  be  sufficient  to  solve  for  the  fracture 
characteristics  (y>  e)  deterministically.  In  that  case,  a  proba¬ 
bility  density  function  p(y,  e)  is  to  be  established  to  identify 
the  probability  of  fractures  at  certain  locations,  such  that 

p(y,  e)  =  p( Y/e)  p(e) .  (5) 

As  an  example,  when  N  =  1  and  M  =  2,  the  solutions  to  frac¬ 
ture  characteristics  form  a  discrete  set  of  K  pairs  such  that 
K  >  1.  The  probability  distribution  function  of  the  pairs  in  the 
solution  set  is  then  given  by: 

K 

F(y,  e)  =  l  p  -H[(  y»  e)  -  G.(R  •)],  i  -  1,  2  (6) 

k=l  J  J  1 

where  H(x)  is  the  Heaviside  function  that  H(x)  =  x  <  0 

Pj  is  the  probability  that  the  solution  equals  the  jth  pair 
(y»  e)j  in  the  solution  set  and 
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(7) 


K 


When  N  *  1  and  M  =  1  there  will  be  infinitely  many  possible  solu¬ 
tion  pairs  to  the  equation  (4),  in  which  case  F(y,  e)  will  be  a 
continuous  function. 
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2.2  Fracture  Hinge 

As  a  point  of  departure,  we  shall  consider  a  fracture 
geometry  shown  in  Figure  la.  Under  flexural  loading,  the  crack 
(a)  is  the  Mode  I  (open-mode)  crack.  The  intensity  of  the  frac¬ 
ture  damage  in  such  a  geometry  is  defined  by  -y  =  a/b  at  the  posi¬ 
tion  e  =  L*/L.  The  existence  of  a  crack  in  a  structural  member 
renders  an  overall  softening  effect  in  the  compliance  of  the 
structure.  The  reduction  in  stiffness  can  be  quantified  by  the 
location  and  the  geometry  of  the  crack.  For  bar  structures,  the 
lower  frequency  modes  are  generally  flexural.  A  crack  in  a  flex¬ 
ural  member  suggests  a  slope-discontinuity  at  the  section  con¬ 
taining  the  crack,  although  the  displacement  is  still  continuous 
there.  The  effect  of  fracture  at  a  section  may  therefore  be 
represented  by  a  spring-loaded  hinge  of  a  torsional  spring  coef¬ 
ficient,  k,  which  may  be  determined  either  analytically  or  exper¬ 
imentally  [5],  Figure  lb.  The  principal  concerns  are  whether  the 
spring  coefficient  is  a  constant  (i.e.,  load/displacement- 
independent)  and  is  invariant  to  the  location  along  the  flexural 
member.  In  other  words,  for  analytical  feasibility,  the  "frac¬ 
ture  hinge"  should  be  entirely  characterized  by  the  intensity  of 
the  fracture  damage  which  is  y  =  a/b  in  the  present  case. 

Development  of  the  simulation  is  to  be  based  on  the  equiva¬ 
lence  of  energy  contents  of  flexural  members,  one  with  a  cracked 
section,  another  with  a  hinge.  If  the  "fracture  hinge"  in  the 
simulated  beam  is  linear,  the  energy  stored  in  the  torsional 
spring  is 


(8) 


where  M  is  the  resisting  moment  at  the  location  of  the  fracture 
hinge.  It  is  observed  that  when  there  is  no  crack  the  spring 
constant  approaches  infinity.  The  softening  effect  of  the  crack 


plane-strain  h=2b 


1— L*— 


torsional  spring  constant,  k 


(b) 


Figure  1.  Fracture  geometry  and  equivalent 
fracture  hinge. 
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results  in  an  increase  of  energy  content  under  load-prescribed 
condition  [6].  The  increase  in  energy  is  provided  by  the  surface 
energy  of  the  crack.  Hence,  assuming  that  the  crack  geometry  is 
stable  for  Mode  I  crack. 


where  g  is  the  fracture  toughness,  A  is  the  area  of  one  side  of 
the  crack  surface,  Kj  is  the  stress  intensity  factor  for  Mode  I 
crack.  The  stress  intensity  factor  is  related  to  the  loading  and 
geometry  as 

Kj  -  /fa  f( y)  =  jv^a  f(£  )  (10) 

where  a  is  the  crack  depth,  b  is  the  half  depth  of  the  beam,  M  is 
the  resisting  moment  at  the  cracked  sction,  Z  is  the  section 
modulus  and  f(Y)  is  the  dimensionless  fracture  intensity  factor. 
For  a  beam  of  rectangular  cross-section,  f(y)  can  be  obtained 
from  expressions  given  by  Bentham  and  Koiter  [7]  as 

f( y)  ■  (1  -  y)"1,5  [1.122  -  2.363y  +  4.367T2  -  4.88  y3 

+  2.845  y4  -  0.663  y5].  (11) 

Combining  equations  (8,  9,  10)  we  thus  arrive  at  the  spring 
constant 

a[f(Y)]2dA.  (12) 


For  a  rectangular  cross  section,  Z  *  2/3  b^w,  where  w  is  the 
width.  Therefore 
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<  = 


(13) 


2Eb2w  /ry  x[f(x)]2dx 

9(1  -  v2)it /Jo 

or  in  its  dimensionless  form  (shown  in  Figures  2a,  b,  c) 

-K.iL^p«,z/r  X[f*)]2d*  (13a) 

EbS  y  Jo 

The  expression  in  (13)  therefore  shows  that  the  torsional  spring 
coefficient  is  entirely  determined  by  the  crack  geometry  for  a 
locally  uniform  beam  of  a  specific  material.  The  equivalent 
"fracture  hinge"  indeed  has  a  linear  torsional  spring. 

It  is  noted  that  Figure  la  illustrates  a  beam  with  symmetri¬ 
cal  cracks.  During  flexural  vibration,  the  tension  side  will 
open  to  provide  the  softening  effect.  The  compression  side  is 
closed,  thus  contributes  nothing  toward  the  hinge  effect.  The 
fracture  hinge,  therefore,  has  a  linear  torsional  spring  provided 
that  the  cracks  on  two  sides  are  symmetrical.  Asymmetry  leads  to 
bilinear  springs.  In  that  case  the  torsional  spring  will  be 
approximated  by  a  linear  spring  with  equivalent  period.  As  an 
illustration,  we  shall  consider  the  extreme  case  of  a  one-sided 
crack.  In  such  case,  the  frequency  of  the  beam  is  the  same  as 
that  of  an  uncracked  one  when  the  crack  is  closed;  the  frequency 
will  be  that  of  a  symmetrically  cracked  one  when  the  crack  is 
open.  Let  t^/2  be  the  half -period  of  modal  oscillation  when 
the  cracked  side  is  in  compression  and  tj/2  be  the  half-period 
when  the  cracked  side  is  in  tension.  Figure  3  shows  the  one 
complete  period  of  a  one-side  cracked  beam.  The  period  of  the 
beam  with  one-side  crack  is  therefore 

t  =  |  (ti  +  t2). 


13 


Equivalent  torsional  spring  constant  k  =  (1  -  v2)ic/Eb2w 


Fracture  damage  y  =  a/b 


Figure  2a.  Spring  constant  for  shallow  crack. 
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Equivalent  torsional  spring  constant  k  =  (1  -  v2)ic/Eb2w 

...  m  01  sianfl-1  _ <-*->  -p.  ui  cn 


Fracture  damage  y  =  a/b 
Figure  2b.  Spring  constant  for  deeper  crack 


Period  of  a  one-sided  crack. 
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When  we  designate  o»i  the  modal  frquency  of  the  uncracked  beam. 
Awl  its  change  for  beam  with  symmetrical  cracks,  and  aoj  the 
frequency  change  of  the  one-side  cracked  beam,  we  have  then 
=  2n/u>p  T2  =  2ir/(u)^  -  Aui^)  and  t  =  2ir/(w^  -  Ao>),  or 


Frequency  variation  IT  is  therefore 

Tf  -  4* - ^ -  or  R  =  — 

(2  -  R)  (1  +  R) 

where  R  is  the  frequency  variation  of  beam  with  symmetrical 
cracks  for  the  same  depth  of  crack.  Figure  4  shows  the  relation¬ 
ship  of  the  frequency  variations.  The  probability  of  asymmetry 
at  the  cracked  section  accounts  for  the  random  characteristics  of 
the  function  in  (4).  In  the  present  report,  the  fracture  diagno¬ 
sis  will  be  limited  to  that  of  symmetrical  cracks  only.  The  pair 
of  fracture  damage  characteristics  are  then  y  =  a/b  and  e  s  L*/L. 
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Frequency  variation  of  one-side  cracked  beam 


0.5  1.0 


R  =  Frequency  variation  of 

symmetrically  cracked  beam 


Figure  4.  Frequency  variations. 


2.3  Fracture  Diagnosis 

In  order  to  identify  the  damage  on  the  structure,  it  is 
important  to  have  a  priori  a  definitive  effect  of  the  damage  on 
the  structure;  that  is,  it  is  necessary  to  establish  the  damage 
function  (4),  analytically  or  experimentally.  For  structures 
composed  of  flexural  members,  the  deformation  of  the  entire 
structure  can  be  defined  by  a  lateral  displacement  set 
{wi(xi,t) },  governed  by  a  set. of  differential  equations 

4  2 

3  w,  3S*. 

(El). - J+  p.  — -!-  0  (15) 

3X^  1  3t^ 

where  (EI)i  is  the  beam  stiffness  and  is  the  lineal  mass 
density  of  the  ith  section  of  the  structure;  summation  convention 
is  not  applied  in  the  equation.  At  the  modal  oscillation  of 
frequency  id,  the  displacement  is  defined  by  the  modal  shapes 
{yi(xi) } 

^(x. ,  t)  =  y^  (x. )  sin  <ot. 

The  set  of  differential  equations  (15)  become  the  modal  shape 
equations; 


to  p  . 


TtrjY 


where  the  primes  (')  indicate  derivatives  with  respect  to 
the  argument  x,.  In  terms  of  dimensionless  variables 

5.  =  x^/L^,  *  u>^p.L^/(EI) .  ,  the  modal  shape  equations 

become 


(16) 
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-  ■  *■ 


Let  there  be  a  reference  section  (p,  El,  L) ;  the  characteristics 
of  all  other  sections  can  be  defined  as 


4 


(0/6) 4 


fi  El 
p  ( E I )  i 


ni  =  (EI)-j/(EI)  and  t-j  =  L/L-j .  The  general  solutions 
are: 


(17) 


y1(C1)  =  Aicosh(aiS5i)  +  Bisinh(aiSCi)  +  Cicos(«135i) 

+  Disin(o1BC1).  (18) 

With  application  of  boundary  conditions,  a  characteristic 
equation  can  be  established 

g/0)  -  0.  (19) 

From  the  set  of  characteristic  values  (0}  the  modal  frequencies 
(u}  can  be  determined  as 


When  a  crack  is  developed  at  a  location  L*,  as  shown  in  Figure 
la,  the  member  is  to  be  divided  into  two  separate  flexural  mem¬ 
bers  with  modal  shapes  {yj,  y£ }  that  satisfy  certain  continu¬ 
ity  conditions  at  the  crack  location.  The  conditions  are 

y/L*)  =  y2(l*),  y/'U*)  *  y/'U*),  yf'd*)  •  y2'"(L*) 


and 


<[y2'(L*)  -  y/a*)]  =  EIy:  ' ' (L*) . 


21 


I 


The  first  three  conditions  are  displacement  and  forces  continu¬ 
ity,  the  last  condition  is  the  spring  hinge  effect.  In  terms  of 
dimensionless  quantities,  e  =  L*/L  and  9  =  EI/kL,  the  conditions 
become 

•N 

y^e)  =  y2(e),  y^  1  (e)  =  y2‘ '  (e),  y^  1  •  (e)  =  y2*  "  (e) 


and 


ey^'Ce)  =  y2'(e)  -  yx* (e)  J 


The  weakened  structure  will  result  in  another  characteristic 
equation 


g2(  01  9»  e)  =  0.  (22) 

where  9  and  e  are  the  governing  parameters  for  the  solution  set 
of  new  characteristic  values,  [0(9,  e)].  The  new  frequency  set 
is  similar  to  (20)  as 


w(0,e) 


a)  -  Aui(  9,e) 


(20a) 


From  the  expressions  in  (20)  and  (20a),  the  damage  functions  (4) 
are  determined.  The  pair  of  governing  parameters  are  9  and  e, 
rather  than  y  and  e.  For  a  rectangular  cross  section,  9  has  the 
value. 


9  =  3ir(l  -  v2) 


(r)/V<A)]  2d* 


(23) 


The  parameter  e  is  seen  to  be  principally  affected  by  the  frac 
ture  damage  (y)  and  the  slenderness  ratio  (b/L).  The  measure¬ 
ments  of  frequency  variations  {R-j }  are  related  to  fracture 
damage  through  this  dimensionless  parameter,  hence  called  the 
sensitivity  number;  that  is,  the  greater  the  value  of  the  sens 
tivity  number,  the  larger  are  the  variations  in  modal  frequen¬ 
cies.  Its  value,  in  e/(l  -  ,  is  graphed  in  Figure  5. 
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3.0  DAMAGE  DIAGNOSIS  IN  SIMPLE  STRUCTURES 

The  principle  of  damage  diagnosis  is  to  be  applied  and 
illustrated  to  some  simple  structures;  namely  a  simple  beam,  an 
L -frame,  and  a  U-frame,  Figure  6. 

3.1  Simple  Beam 

The  modal  frequencies  for  a  simple  beam  can  be  readily 
established  analytically.  Those  of  an  uncracked  beam  are  avail¬ 
able  from  many  standard  publications,  such  as  Young  and  Felgar 
[8].  The  solution  of  the  shape  function  y(s)  is 

yU)  =  A  cosh(sc)  +  8  sinh(ec)  +  C  cos(ec)  +  D  sin(eO  (24) 


Two  homogeneous  kinematic  forced  boundary  conditions  or  two 
kinetic  natural  boundary  conditions  or  mixed  two  are  available  at 
each  end  (c  =  0  o r  1).  For  the  resulting  set  of  homogeneous 
equations,  the  determinant  of  coefficients  for  (A  8  C  D)  must  be 
zero  for  non-trivial  solutions.  The  characteristic  equation  in  e 
(19)  is 


g^B)  » 


U6(B) 


(19a) 


where  lli(A)  and  Ug(B)  are  two  (2  x  4)  matrices  corresponding 

to  the  boundary  conditions  at  A  and  B,  i.e.,  £  *  0  and  1,  respec¬ 
tively.  For  instance,  for  the  cantilever  end  conditions 


yx(A)  « 


1 

0 

L. 


0 

1 


1  0 
0  1 


(25) 


U6(B)  = 


coshe  sinhe 
sinhe  coshe 


-cose 

sine 


-sine 

-cose 


(26) 
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The  solution  set  {g^}  thus  leads  to  the  modal  frequencies 
{a>i }  of  an  uncracked  beam  by  (20). 


For  a  beam  with  a  crack  at,  the  section  L*,  Figure  la,  the  modal 
shape  function  is  divided  into  two  parts,  yiU)  and  y2(0, 
for  the  range  (0,  e)  and  (e,  1),  respectively,  where  e  =  L*/L. 
Each  portion  of  the  shape  function  must  satisfy  the  governing 
equation  (16)  and  have  a  general  expression  as  (18).  The  addi¬ 
tional  continuity  conditions  (21)  expand  the  order  of  determinant 
of  the  characteristic  equation  by  four,  to  result  in 


g2(e)  3 


U_1(A)  9 
u3(e)  u4(o) 

0  U6(B) 


=  0 


(27a) 


where  both  null  matrices  0  are  (2x4),  113(e)  and  114(e)  are 
both  (4  x  4)  with  expressions: 


y3(e,e,i) 


0  0  0 

0  0  0 

0  0  0 

96  cosh(lie)  eg  sinhCee)  -9?  cos(ie) 


0 

0 

0 

-es”  sin(?e) 


-  U4(6,e) 


(27) 
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U4(e,e)  = 


cosh  (Je)  sinh  (Je)  cos  (Je)  sin  (Je) 

cosh  (Je)  sinh  (Je)  -cos  (Je)  -sin  (Je) 

sinh  (Je)  cosh  (Je)  sin  (Je)  -cos  (Je) 

sin  (Je)  cosh  (Je)  -sin  (Je)  cos  (Je) 


(28) 


The  frequency  variations  {R^}  can  be  computed  as  functions  of 
(0,  e) ,  subsequently  (y,  e;  b/L),  where  the  slenderness  ratio  b/L 
is  an  influencing  parameter.  Figure  7  (a,  b,  c)  show  three  modal 
frequency  variations  for  a  cantilever  beam  of  a  slenderness  ratio 
of  b/L  =  0.05,  representing  the  function  f,  in  (3).  Figure  8 
shows  the  corresponding  inverse  transform  for  the  damage  function 
G|(  in  (4),  with  specific  given  frequency  variations: 

Rj  3  0.02,  R2  =  0.01,  and  R3  3  0.0225.  For  a  beam  with  a 
single  crack,  the  two  unknowns  (y,  e)  should  be  solvable  with  two 
frequency  variations  (say  Rj,  R2) .  However,  as  it  is 
observed  in  Figure  8  for  the  range  of  e:  [0,  1],  there  are  two 
solution  points:  e][  *  0.08  and  e2  *  0.33,  the  intersections 
of  Ri  and  R2  curves.  Yet,  the  solution  must  also  satisfy  the 
other  frequency  variations  as  well;  in  the  present  case,  the 
third  mode.  The  only  solution  is  therefore  at  e  =  0.33  corres¬ 
ponding  to  a  crack  damage  intensity  of  0.32.  Figure  9  shows  the 
damage  function  (4)  for  given  modal  frequency  variations: 

Rl  3  0.026,  R2  *  0.096,  and  R3  3  7.33  x  10'5.  The 
resulting  fracture  characteristic  pair  is  (y,  e)  *  (0.6,  0.5). 


Third  node  frequency  yertetUm  (Rj) 


Crtc*  position  (•  ■  L*/l) 

Flguro  9.  Outgo  function  (r.o)  •  S(R,) 


3.2  L-Frame 

The  modal  frequency  analysis  of  a  simple  L-frame  is 
usually  carried  out  as  two  beams  rigidly  connected  at  the  cor¬ 
ner.  The  contraction  in  axial  length  is  treated  as  second-order 
smallness  in  comparison  to  the  lateral  displacement.  In  the  case 
when  both  anchors  are  either  fixed  or  pinned,  the  displacement  at 
the  junction  corner  is  therefore  approximated  by  order  magnitude 
to  be  zero.  If  one  of  the  anchors  is  allowed  to  displace  axi¬ 
ally,  the  inertial  effect  of  the  beam  containing  the  free  termi¬ 
nus  must  be  accounted  for  by  Newton's  second  law.  Since  the  free 
terminus  condition  will  be  demonstrated  in  the  U-frame  problem, 
we  shall,  for  illustrative  purpose  of  the  damage  diagnosis, 
choose  as  an  example  an  L-frame  with  both  termini  fixed.  Without 
loss  of  generality,  the  cracked  section  is  placed  on  the  vertical 
leg.  Figure  10.  The  continuity  conditions  and  geometrical  con¬ 
straints  at  the  corner  with  reference  to  the  general  solution 
(18)  are: 

yi(l)  •  0,  y2  (0)  *  0,  iiyi'  (1)  =  i2  y2'(0),  (29) 

*1  Vl',(1)  =  *2  n2^2* 

where  =  L/L-j ;  and  n-j  =  (EI)-j/(EI) .  The  characteristic 
equations  (19  and  22)  of  the  uncracked  and  the  corresponding 
cracked  l-frames  are: 

o 

U5(B)  -  0  (19b) 

i(« 

ux(A)  o  o 

U3(e)  U4(e)  0 

92(^  3  o  u2(b)  u5(b) 

o  "  o  u6(C) 


34 


y2 


Figure  10.  L-Frame  with  a  single 
cracked  section. 
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where  0  and  J  are  the  characteristic  values  of  the  uncracked  and  the 
corresponding  cracked  L-frames,  such  that  s4  =  w2pL4/(EI)  as 
referred  to  the  reference  beam,  the  vertical  leg.  Ui(A)  and  Ug(C) 
are  (2  x  4)  matrices,  related  to  the  anchor  conditions  at  A  and  C; 
specifically  for  fixed  ends  as  shown  in  Figure  10,  Ui(A)  has  thus  the 
same  expression  as  (25),  Ug( C)  differs  from  (26)  since  the  end 
conditions  are  different,  that  is 


U_i(A)  = 

'101 

0 

U  ■) 

010 

1  . 

u6(C)  = 

'  cosh(<x20) 

sinh((*20)  cos(<x20) 

sin((*20) ' 

(30) 

-u 

sinh^e) 

cosh(ot2&)  -sin(a20) 

COS(ct20) 

where  ai-  =  0 \/ Q.  ^(B)  and  (15(8)  are  (4  x  4)  matrices 

resulting  from  the  continuity  conditions  and  constraints  at  the 
corner  (2°),  that 


U2(B) 


U7( e;l,l)  r-U7(0;O.2) 

Ug(0;l,l)  and  Ug(B)  =  0 

0  Ug(0;O,2) 


where 

I"  a,£,[sinh(a,05)  cosh(a,0?)  -sin(a.0£)  cos (a?0S)  (31) 

u,{  b;c.D  3  I  1  i  2  ,  „ 

I  n.a^e‘[cosh(a.0C)  sinh(ai0£)  -cos(ai0£)  -sinfo^O] 


and 


Ug(0;C,i)  =  [cosh(ai0?)  sinh(ai0C)  cosfa^e)  sin(a.0c)] 


(32) 


where  =  (EI)i/(EI)  and  £.  =  L/L- .  ( e)  and  U4(e)  are 

unchanged  as  in  the  simple  beam,  (27,  28).  The  solutions  of 
equations  (19b,  22b)  are  substituted  into  (20)  to  solve  for  the 
modal  frequencies  of  the  uncracked  and  the  corresponding  cracked 
L-frames,  {u-j }  and  (w-j  -  Acu-j }  respectively.  Figure  11 
illustrates  the  measurements  of  Rj  =  0.0363  and  Rj  =  0.031 
for  an  L-frame  with  *2  =  1-333,  n2  =  1,  and  <*2  3  0.75, 
using  the  vertical  leg  as  the. reference  member,  and  a  slenderness 
ratio  of  0.05.  The  solutions  of  e  are  multivalued  at  .140,  .155, 
.480,  and  .5jG.  If  R3  =  0.0109  is  also  available,  the  pair  of 
damage  characteristics  are  then  determined  at  y  =  a/b  =  0.4  and 
e  =  .480.  It  is  noted  that  the  function  g2(?)  will  be  changed 
when  the  cracked  section  is  in  the  horizontal  leg,  instead  of  the 
vertical  leg.  The  same  expression  as  in  (22b)  results  when  the 
designations  of  the  horizontal  and  vertical  legs  are  reversed, 
resulting  in  a  new  03(0*)  *  0.  In  other  words,  the  set  of  [0} 
corresponds  to  (y,  e)  on  the  vertical  leg;  and  the  set  of  {e*} 
corresponds  to  (y*»  e*)  on  the  horizontal  leg.  A  given  measure 
of  {Ri  },  therefore,  will  be  consistent  with  only  one  pair  of 
the  damage  characteristics.  Hence,  the  frequency  change  measure¬ 
ment  determines  the  location  (e)  and  the  intensity  (y)  of  the 
damage  uniquely  in  an  L-frame. 

The  unique  solution  will  not  be  available  for  a  special  case 
when  the  expressions  of  g2(7)  and  93(0*),  duplicate.  Such  is 
the  case  for  L -frame  of  equal  legs.  It  is  shown  in  Figure  12 
that  the  given  measurements  Ri  *  0.0363,  R2  =  0.0061  and 
R3  =  0.046  for  an  L-frame  of  slenderness  ratio  0.05  identify 
the  damage  location  e  =  0.05  and  the  damage  intensity  y  -  0.6  in 
either  one  of  the  legs. 
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Crack  position  e  =  L*/L 


Figure  12 


.  Damage  diagnosis  in  equal  leg  L-frame. 


39 


3.3  U-Frame 


The  modal  analysis  of  a  U-frame  is  similar  to  that  of 
an  L-frame  with  an  additional  beam  section.  Figure  13  shows  the 
geometry  and  possible  crack  location.  In  addition  to  the  end 
conditions  at  two  anchors,  there  are  continuity  conditions  at  the 
corners  and  the  kinetic  equation  defining  the  motion  of  the 
horizontal  beam  which  is  not  constrained  in  axial  motion.  The 
end  conditions  will  be  governed  by  how  the  ends  C  and  D  are 
anchored.  The  corner  conditions  are  the  continuity  conditions 
and  the  geometrical  constraint.  At  the  corner  A  in  the 
dimensionless  expressions: 

Vl'f1)  =  VlV*1)  =  n2*22y2"^’  and 

y2(0)  =  0  (33) 

At  the  corner  8: 

*2y2,(1)  a  *3y3 ' (°) ’  n2A22y2"(l)  =  n3*32y3"(0),  and 
y2(l)  ■  0.  (34) 

The  motions  of  the  horizontal  beam  (AB)  are 

yj(A)  =  -y3(B)  or  y1(l)  »  -y3(0)  (35) 

nl4l3yl",(A)  +  =  -»24n2A23yl(A)  (36) 

or 

nlJil3yr  "(1)  +  a2404n2je23yl^1^  +  n3^3 3y3  "  '  (0)  =  0 

where  I i  *  L/Li ,  n^  *  (El). /(El)  and  e24  =  a>2p2L24/(EI)2.  The 
condition  (35)  is  evident  with  reference  to  the  coordinates  shown 
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in  Figure  13.  To  arrive  at  (36),  we  shall  refer  to  Figure  14  for  cross 
shears,  which  are  the  axial  driving  forces  for  the  horizontal  beam 
(AB).  The  axial  acceleration  of  the  horizontal  beam  shall  be  the  same 
as  the  lateral  acceleration  of  the  vertical  leg  at  either  of  the 
corners  A  and  B;  that  is, 

Vj(L)  +  V3(0)  =  (p2L2)[-yj(L)]. 

At  modal  oscillation,  y^(L)  =  -  oi  y^(L).  Hence  for  V  =  -  Ely11', 
-(EI)1y1",(L)  -  (EI)3y3,',(0)  =  p^^y^L) , 


or 


"^"'(L)  +  0^3* '  '(0)  +  C«2p2L|/(EI)2]C(EI)2/El3yi(L)L23  =  0 

When  expressed  in  dimensionless  form,  the  expression  in  (36) 
results.  For  an  uncracked  U-frame,  there  are  two  pairs  of  end 
conditions  for  the  two  anchors,  and  eight  corner  conditions  at  A 
and  B  for  the  three  sets  of  coefficients  {A,  B,  C,  D}  of  beam 
equation  (18)  and  the  characteristic  parameter  e.  The  resulting 
characteristic  equation  (19)  for  the  uncracked  U-frame  is 


arranged  as 

^(C) 

U2(A) 

0 

U5(A) 

0 

0 

gjU)  = 

Ug(A) 

0 

u7(B) 

0 

Ug(B) 

y9(o) 

=  0 

(19c) 

where  all  submatrices  are  referred  to  the  positions  on  the 
U-frame  for  convenience  of  discussion.  For  instance,  as  C  and  D 
are  the  two  ends  of  the  U-frame,  the  matrices  U^(C)  and  Ug  (D) 
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thus  define  the  anchor  conditions  of  the  frame.  For  fixed  endsas 
shown  in  Figure  13 


yi(c) 


ri  o  i 

0  1  0 

L— 


Ug(°)  » 


”cosh(a30) 

^sinhCo^e) 


:] 

sinh(a3e)  cosCa^e)  sinCa^e) 
cosh(agb)  -sin(a3$)  cos(aj0)  . 


112(A)  and  U5( A)  are  both  (3  x  4)  matrices  which  result  from 
the  continuity  conditions  at  A,  (33),  that. 


U2(A)  « 


-10^ 


and  Ug(A) 


-U10(@;0,2) 

Ujj(0  ',0,2) 


(37) 


where 


u10(s;c5i) 


aiii[sinh(ai05)  cosh(Oi0C)  -sin(<xi ee )  cosfa^e)] 


2.2r 


niai4i[cosh(ai0C)  sinh(Oi0£)  -cosfajK)  -sin(Oi0?)] 

(38r 


and 

0^(855, i)  *  [cosh(a^ 05)  sinh(a^05)  cos(ai0€)  sin(Oi0?)]  (39) 


The  matrices  Ug(A),  U7(B)  and  Ug(B)  are  all  (5  x  4)  as  groupings 

of  the  corner  conditions  at  B,  (34),  and  the  equations  of  motion 
of  the  horizontal  beam  (AB),  (35,  36);  they  are 
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u6(a)  = 


yn(e;i,i) 

u12( e;l,l)Q+  a2*6ry23  Uu(b;1,1) 


Uy(  b)  a 


UioCeU.z) 

Uu(e;l,2) 


Uo(B)  = 


where 


yn(8;0,3) 

u12(b;0,3) 

■u10(b;0,3) 


U12(B;5»i) 


[sinh(a.BC)  cosh^BC)  sinU^C)  -cos^B?)] 

(43) 


Depending  on  the  choice  of  the  reference  member  in  the  frame,  the 
geometry/material  ratios  can  be  simplified.  For  instance,  if  the 
horizontal  member  is  chosen  as  the  reference,  a2  =  *2  =  n2  =  1 
and  B2  *  3-  The  characteristic  equation,  gg( 0^)  *  0,  is 
obtained  after  expanding  the  determinant  (19c)  by  an  order  of  4, 
with  the  introduction  of  the  fracture  hinge  conditions  (21),  or 
the  submatrices  113(e)  and  U^e),  (27,  28).  For  the  case  when  the 
crack  occurs  on  the  vertical  leg,  AC,  we  have 


Ujtc)  0 

u>  y4(e) 

92(b)  *  0  U2(A) 

?  U>> 

0  0 


Ug(A)  0 
U?(B)  u8(B) 

0  Ug(D) 
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When  the  crack  occurs  on  the  horizontal  member,  AB,  the  hinge  is 
between  A  and  B,  resulting  in: 


yi(c)  ooo 

U2(A)  U5(A)  0  0 

g2(0)  =  *  0  U3(e)  U4(e)  0  =  0  (22d) 

U6(A)  "  0  U?(B)  Ua(B) 

"  o  o  "  o  u9(d) 

Notice  that  the  fracture  hinge  matrices  (27,  28)  are  inserted 
where  they  separate  the  submatrices  in  a  colunn  which  are  defined 
at  the  two  end  points  of  the  cracked  member.  To  obtain  the  char¬ 
acteristic  equation,  when  a  crack  is  in  the  other  vertical  leg, 
one  may  either  insert  the  fracture  hinge  matrices  separating 
Ug(B)  and  Ug(D),  or  use  the  similar  expression  of  (22c)  with  a 
left/right  transfer  of  the  configuration.  In  the  expressions 
(22c)  and  (22d),  the  sensitivity  number  I  is  defined  as 

0  =  [EI/kL]1  (44) 

where  the  index  ' i 1  denotes  the  specific  cracked  member  in  the 
structure. 

The  characteristic  equations  (19c,  22c,  and  22d)  are  quite 
complex.  The  solutions  of  (e,  e)  then  (y,  e)  =  G(Ri),  equation 
(4),  are  thus  programmed  in  a  computer  without  curve  plots  such 
as  those  in  the  simple  beams  (Figures  7  through  9)  and  the 
L-frames  (Figures  11,  12).  Both  characteristic  equations,  (22c 
and  22d)  are  used  in  the  diagnosis.  Again  three  measurements  of 
frequency  variation,  Rj,  Rg,  and  R3,  are  sufficient  for  the 
determination  of  the  fracture  characteristics  of  a  single  crack. 
Except,  of  course,  a  U-frame  with  equal  vertical  legs,  which  are 
anchored  the  same  way,  will  have  symmetrical  crack  locations  and 
damage  intensities  that  are  not  differentiable.  In  the  computer 
code,  a  crack  is  first  sought  on  the  vertical  legs  and  then  on 


the  horizontal  member.  In  the  deterministic  case  for  a  single 
crack,  the  three  measurements  of  frequency  variations  will  locate 
the  crack  specifically  on  the  member  of  the  U-frame  where  the 
crack  occurs.  However,  when  inadequate  measurements  are  given, 
probable  fracture  locations  may  occur  on  both  of  the  legs  and  the 
horizontal  member.  For  instance,  as  indicated  by  equations  (6, 
7),  we  have  two  measurements  of  frequency  variations, 

{Rp  R2}  =  {0.022,  0.005} 

for  a  symmetrical  U-frame.  Referring  to  the  horizontl  member  as 
the  reference  beam  we  have 

n2  =  l2  ~  ®2  =  nl  =  n3  =  1*  *1  =  *3  =  1,25>  “1  =  “3  =  0,8 
The  admissible  solututions  are 

{( e,e) }x  ■  {(0.031,  0.050)  (0.426,  0.390} 

{( 9,e) }~  *  {(0.095,  0.078)  (0.180,  0.187)  (0.180,  0.813) 
c  (0.095,  0.922)} 

{( 9,e) }3  -  {(0.426,  0.610)  (0.031,  0.950} 

Because  of  the  symmetry  condition,  there  are  only  four  pairs  of 
solutions.  They  are 

{(0,e)}1  *  {(0.031,  0.050)  (0.426,  0.390)} 

{(e,e)}2  =  {0.095,  0.078)  (0.180,  0.187)} 

as  shown  in  Figure  15.  With  reference  to  equations  (6,  7),  it  is 
proposed  that  the  crack  occurs  at  the  most  likely  location,  e, 
where  the  largest  probability  function,  pj,  occurs.  It  is 
assumed  that  the  largest  pj  of  the  set  in  turn  corresponds  to 
the  location,  e,  where  the  moment  distribution  has  its  greatest 
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Symmetry  plane 


(0.180,  0.187) 
(0.095,  0.078) 

r — » — 


\ 

I 


(0.426,  0.390) 


(0.031,  0.050) 

/W 


Ri  =  0.022  Rz  =  .005 


Figure  15.  Fracture  damage  characteristic  pair  (Q,e) 
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value  among  the  possible  solutions.  Figures  16  and  17  illustrate 
the  modal  shapes  and  the  corresponding  moment  distributions  of 
the  symmetrical  l)-frame.  For  illustrative  purposes  we  may  assume 
that 

:n,j  • 

The  approximation  indicates  that  in  the  present  example,  since 

Rj  is  more  than  four  times  R2,  the  resisting  moment  M(ej) 

is  to  be  based  only  on  the  first  mode  moment  distribution.  Hence 

{Pj}  =  {mj}  =  {0.45  0.15  0.17  0.23} 

The  most  probable  location  is  therefore  on  either  of  the  vertical 
legs  at  a  distance  of  0.05  Li  from  the  anchor.  The  fracture 
damage  intensity,  y,  corresponding  to  the  sensitivity  number, 
e  a  0.031,  depends  on  the  slenderness  ratio  (b/l)?  re:  equation 
(23).  If  we  use  a  rectangular  cross  section,  for  a  5  percent 
slenderness  ratio,  the  damage  intensity  is  y  =  0.33;  for  a  1.5 
percent  slenderness  ratio,  it  takes  a  damage  intensity  of 
y  *  0.52  to  effect  a  corresponding  frequency  change.  The  values 
of  y  can  be  approximately  deduced  by  referring  to  Figure  5. 
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moment  distribution  on 
vertical  elements 


— — —  —  on  the  horizontal  element 


Figure  17.  Moment  distributions  at  modal  oscillation 
in  a  symmetrical  U- frame. 

-  £3  1.25,  n x  =  n 3  =  1,  ctx  =  ct$  =  0.8 
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4.0  CONCLUSION 

The  present  report  addresses  the  effect  of  fracture  damage 
on  the  stiffness,  then  on  the  modal  frequencies  of  a  structure. 
The  intensity  of  damage  and  its  location,  (y,  e),  can  be  identi¬ 
fied  if  sufficient  measurements  of  frequency  variations,  {R-j}, 
and  the  damage  function,  Gp(R-j),  can  be  established.  For 
simple  structures  such  as  uniform  beams,  L-frames  and  U-frames, 
the  analytical  form  of  the  damage  function  can  be  readily  estab¬ 
lished  or  programmed  in  computers.  The  deterministic  solution 
depends  on  the  adequate  number  of  measurements  in  frequency  vari¬ 
ations.  When  an  insufficient  number  of  measurements  is  avail¬ 
able,  the  diagnosis  can  still  be  carried  out  on  a  probabilistic 
base.  At  present,  however,  there  is  no  quantitative  measure  of 
probability,  since  the  postulation  of  using  moment  distribution 
as  an  analogy  to  the  probability  distribution  is  at  best 
qual itative. 

The  introduction  of  the  fracture  hinge  indeed  simplifies  the 
mathematics  of  damage  function.  The  assumptions  and  the 
subsequent  development  are  based  on  the  concept  that  the  strain 
energy  increment  under  constant  load  comes  from  the  total  conver¬ 
sion  from  the  surface  energy  that  results  from  creation  of  the 
crack  surface.  The  very  simplicity  is  the  beauty  of  the  theory. 
The  extent  of  deviation  from  the  physical  behavior  of  the  struc¬ 
ture  must  yet  be  established  by  experiment.  Further  work  may 
then  be  expanded  to  cracks  of  different  geometry. 

It  is  noticed  that,  while  the  fracture  damage  is  character¬ 
ized  by  the  pair  (y,  e) ,  y  *  a/b  and  e  =  L*/l,  the  damage  inten¬ 
sity  is  really  measured  by  the  sensitivity  nunber,  e.  Indeed, 
the  greater  is  the  fracture  damage  intensity,  the  greater  is  the 
sensitivity  number;  but  for  a  given  y,  the  more  slender  the  beam, 
the  smaller  is  e.  From  a  practical  point  of  view,  the  sensitiv¬ 
ity  number  is  indicative  of  how  easy  or  difficult  it  is  to  detect 
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a  crack.  A  very  slender  beam  may  have  a  severe  crack  on  it 
( i * e - »  y  large),  yet  changes  in  the  modal  frequencies  would  be 
small,  thus  difficult  to  detect. 

In  a  fractured  structure,  the  modal  frequency  variations 
are,  of  course,  influenced  by  the  location  of  the  crack, 
e  =  L*/L.  The  sensitivity  measure  should  seem  to  be  a  number 
combining  both  0  and  e.  However,  the  fracture  hinge  matrices 
^3(3*  e,  9)  and  U^(b;  e) ,  equations  (27,  28)  in  the  frequency 
equation  (22),  demonstrate  that  the  two  quantities  must  be 
treated  as  two  separate  entities.  It  is  this  very  property  that 
allows  determination  of  both  the  location  and  intensity  of  the 
fracture  damage. 


REFERENCES 


1.  Liebowitz,  H.,  "Fracture  and  Carrying  Capacity  of  Notched 
Columns,"  Chap.  4,  Fracture,  Vol.  IV,  ed.,  Liebowitz,  Acad. 

Press,  1969,  pp.  113-171. 

2.  Bennett,  J.  G. ,  Ju,  F.  0.,  and  Anderson,  C.  A.,  "An  Investi¬ 
gation  of  Failure  Mechanisms  for  HTGR  (High  Temperature 
Gas-cooled  Reactor)  Core  Supports,"  LASL  Report 
LA-NUREG-6639-MS  (Dec.  1976). 

3.  Ju,  F.  D.,  Bennett,  J.  G.,  and  Anderson,  C.  A.,  "Structural 
Safety  Analysis  of  HTGR  Core  Supports,"  Advances  in  Civil 
Engineering  through  Engineering  Mechanics,  Proc.  2nd  EMD 
Specialty  Conference,  ASCE,  1977,  pp.  427-430 . 

4.  ASME  B  and  PV  Code,  Sec.  XI,  ART.  IWA-2000,  1974. 

5.  Chondros,  I.  G. ,  and  Dimarogonas,  A.  D.,  "Identification  of 
Cracks  in  Welded  Joints  of  Complex  Structures,"  Jour.  Sound 
Vibra.,  1980,  pp.  531-538. 

6.  Tada,  H.,  Paris,  P.,  and  Irwin,  G.,  "Compliance  Calibration 

Methods,"  The  Stress  Analysis  of  Cracks  Handbook,  Del  Res. 
Corp.,  1973;  pp7  l’.T0-O2  - ' 

7.  Bentham,  J.  P.,  and  Koiter,  W.  T.,  "Asymptotic  Approximation 

to  Crack  Problems,"  Methods  of  Analysis  and  Solutions  of 
Crack  Problems,  Mechanics  of  Fracture  1,  ed.,  G.  C.  Sih, 
r973rpp:"T55-T58. - 

8.  Young,  D.,  and  Felgar,  R.  P.,  Jr.,  Tables  of  Characteristic 
Functions  Representing  Normal  Modes  of  Vibration  of  a  Beam, 
Eng.  Res.  Bui.  No.  4913,  BER,  University  of  Texas,  July 
1949. 


54 


